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We report the observation of TeV 7-rays from the Cygnus region using the 
ARGO-YBJ data collected from 2007 November to 2011 August. Several TeV 
sources are located in this region including the two bright extended MGRO 
J2019+37 and MGRO J2031+41. According to the Milagro data set, at 20 TeV 
MGRO J2019+37 is the most significant source apart from the Crab Nebula. No 
signal from MGRO J2019+37 is detected by the ARGO-YBJ experiment, and 
the derived flux upper limits at 90% confidence level for all the events above 600 
GeV with medium energy of 3 TeV are lower than the Milagro flux, implying that 
the source might be variable and hard to be identified as a pulsar wind nebula. 
The only statistically significant (6.4 standard deviations) 7-ray signal is found 
from MGRO J2031+41, with a flux consistent with the measurement by Milagro. 

Subject headings: gamma rays: general — pulsars: individual (MGRO J2019+37, 
MGRO J2031+41) 



1. Introduction 



The Cygnus region is t he brightest diffus e 7-ray emitti ng region in the no rthern sky 
as observed by both Fermi ( lAbdo et al.ll201ll ) and EGRET (IHunter et al.lll997l ). Complex 
features have been observed in the wavelength bands of radio, infrared, X-rays, and 7-rays. 
This region is rich in potential cosmic-ray acceleration sites, e.g., Wolf-Rayet stars, OB 
associations and supernova remnants. Recently, 24 7-ray sources, including 7 pulsars, have 
been detected using Fermi Large Area Telescope (L AT) two-year data within the region with 
65° < I < 85° and —3° < b < 3° (lAbdo et al.ll201ll ). These are considered candidate sources 
of very high energy (VHE) 7-rays. The Cygnus region is, therefore, a natural laboratory to 
study the origin of cosmic rays. 

Several VHE 7-ray sources have been detected within the Cygnus region in the past 
decade. Th e first was TeV J2032+4130, discovered by th e HEGRA collaboration (lAharonian et al 



20021. 



20051) and con firmed by the experiments Whipple ( Konopelko et al. 2007 ) and MAGIC 



([Albert et al.ll2008l ). Its extension is estimated to be about 0.1°. The power-law spectral 
index is about —2.0 and the integral flux above 1 TeV is 3—5 % that of the Crab flux. MGRO 



J2031+41, detected by the Milagro experiment at 20 TeV flAbdo et al.ll2007al ). is spatially 
consistent with the source TeV J2032+4130, while the measured extension is much larger 
with a diameter of 3.0°±0.9°. This source is likely to be associated with the pulsar 2FGL 



90146 Palermo, Italy. 
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J2032.2+4126, detected by Fermi flAbdo et all 120111 ) 



Evidence of TeV emission at the 4.0 standard deviation (s.d.) level after a n X-ray flare 



from Cyg X-l was observed by the MAGIC experiment on 2006 September 24 (jAlbert et al. 
20061 ). 



The source VER J 2019+407 was d iscovered in a survey of the Cygnus region by the 
VERITAS experiment (IWeinstein 1120091 ). The measured extensions are 0.16° ± 0.028° and 
0.11°±0.027° along the major and minor axes, respectively. This source is spatially coincident 
with the Fermi source 2FGL J2019.1+404 0, which is potenti ally associated with a supernova 
remnant or a pulsar wind nebula (PWN) (jAbdo et al.ll201ll ). 



During the deep VERITAS observations of the C yg OBI regi on, a point source VER 
J2016+372 was discovered at the location of CTB 87 (kliu I boilh . The flux is about 1% 
that of the Crab and the spectral index is about —2.1 from a preliminary analysis. 

This region also contains the bright unid entified source MG RO J2019+37, which was 



detected by the Milagro experiment at 20 TeV (jAbdo et al.ll2007al ) and is the most significant 
source in the Milagro data set apart from t he Crab Nebula. It s extension is a = 0.32°±0.12° 



in a symmetric two-dimensional Gaussian (jAbdo et al.ll2007b). which has 68% o f the events 



contained in a region with an angular diameter of 1.1°±0.5° (jAbdo et al.ll2007al ). The spec- 



trum of 


this s 


( Smith 


2010) 



At the location of MGRO J2019+37, a 2. 2 s.d. signal correspon ding to 



30% Crab unit was observed by the Tibet AS7 exper iment (lAmenomori et al.ll2010l ). How 



ever, about 0.9° away, a possible source was detected (jAmenomori et al.l 120071 ). The MGRO 



J2019+37 is spatially coincident with the Fermi source 2FGL J2018. 0+3626 and pulsar 



2FGL J2021.0+3651 (jAbdo et all 120111 ) . VERIT AS has surveyed this region, but no emis- 



sion from MGRO J2019+37 has been detected (jWeinstein 1 120091 ). Recently, an in-depth 



observation of the Cyg OBI region has been c arried out b y VERITAS, which unveiled com- 



plex TeV emission around MGRO J2019+37 (lAliu 1 120111 ). 



Among the four known VHE 7-ray sources inside the Cygnus region, MGRO J2019+37 
is enigmatic due to its high flux not confirmed by other VHE 7-ray detectors. The mea- 
surement of the energy spectrum or an upper limit around several TeV is therefore very 
useful to understand the nature of the source and its emission mechanism. The ARGO-YBJ 
experiment is an air shower array with large field of view and can continuously monitor the 
northern sky. The total exposure to the Crab Nebula reaches about 1200 days and its photon 
flux has been detected with a statistical significance of 17 s.d. at energies around 1 TeV, 
which is com parable with the e ight-year value of 17.2 s.d. obtained at energies around 35 TeV 
by Milagro (jAbdo et al. 2009 ). This work presents the observation results for the Cygnus 
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region, including sources MGRO J2031+41 and MGRO J2019+37, with the ARGO-YBJ 
experiment. 



2. The ARGO-YBJ experiment 



The ARGO-YBJ experiment, located in Tibet, China, at an altitude of 4300 m a.s.L, 
is the result of a collaboration among Chinese and Italian institutions and is designed for 
VHE 7-ray astronomy and cosmic ray observations. The detector consists of a single layer of 
Resistive Plate Chambers (RPCs), which are equipped with pick-up strips (6.75 cm x 61.80 
cm each). The logical OR of the signal from eight neighboring strips constitutes a logical 
pixel (called a "pad" ) for triggering and timing purposes. One hundred thirty clusters (each 
composed of 12 RPCs) are installed to form a carpet of about 5600 m 2 with an active area of 
~93%. This central carpet is surrounded by 23 additional clusters ("guard ring") to improve 
the reconstruction of the shower core location. The total area of the array is 1 10 m x 100 



m. More det ails about the detector and the RPC performance can be found in ( lAielli et al. 
20061 . l2009bh . 



The ARGO-YBJ detector is operated by requiring the number of fired pads (Npad) at 
least 20 within 420 ns on the entire carpet detector. The high granularity of the apparatus 
permits a detailed space-time reconstruction of the shower profile and therefore of the inci- 
dent direction of the prim ary particle. The a rrival time of the particles is measured with a 
resolution of about 1.8 ns (jAielli et al.ll2009bl ). In order to calibrate the 18,360 TDC chan- 



nels, we have developed a method using cosmic ray showers flHe et al.ll2007). The calibration 



precision is 0.4 ns and the procedure is applied every month (jAielli et al. 



2009af ). 



The central 130 clusters began taking data in 2006 July, and the "guard ring" was 
merged into the DAQ stream in 2007 November. The trigger rate is 3.5 kHz with a dead 
time of 4% and the average duty cycle is higher than 86%. The angular resolution, pointing 
accuracy and stability of the ARGO-YBJ detector array have been th oroughly tested by 
measuring the shadow of the Moon in cosmic rays (IBartoli et al.ll2011bl ). The point-spread 
function (PSF) is quantified using a parameter i/j 7 q as the opening angle contai ning 71.5% of 
the ev ents. For N pa d > 1000, ip70 is 0.47°, while at N pa d ~ 20 ipm becomes 2.8° (IBartoli et al. 



2011al lbh. This measured angular resolution refers to cosmic ray induced air showers, and it 
is smaller by 30—40% for 7-rays. 
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3. Data analysis 



The ARGO-YBJ data used in this analysis were collected from 2007 November to 2011 
August. The total effective observati on time is 1182.0 d ays. To achieve a good angular 
resolution, the event selections used in ( IBartoli et al.ll2011af ) are applied here and only events 
with zenith angle less than 50° are used. The total number of events after filtering used in 
this work is 1.97xlO n . The opening angle i/j^o for events with N pa d > 60 is 1.36°. In order 
to obtain a sky map, an area centered at the source location in celestial coordinates (right 
ascension and declination) is divided into a grid of 0.1° x 0.1° bins and filled with detected 
events according to their reconstructed arriva l directions. To extr act the excess of 7-rays 
from each bin, the "direct integral method" (IFleysher et al.l 120041 ) is adopted in order to 
estimate the number of cosmic ray background events in the bin. To remove the effect of 
cosmic ray anisotrop y on a spatial scale of 11° x 11°, a correction procedure described in 
( IBartoli et al.l 1201 laf ) has been applied. To take into account the PSF of the ARGO-YBJ 
detector, the events in a circular area centered on the bin with an angular radius of L3j^zq, 
are su mmed together with a weight of the Gaussian-shaped PSF. Equation (17) in (ILi &: Ma 
19831 ) is used to estimate the significance of the excess in each bin. 



With this data analysis, the significance of the excess observed from the direction of the 
Crab Nebula is 17 s.d., which indicates tha t the cumulative 5 s.d. sensitivity of ARGO-YBJ 
has reached 0.3 Crab unit for point sources ( jCao fc Chenll201ll ). For an extended source with 
a symmetric two-dimensional Gaussian shape with a = 0.32°, the sensitivity is degraded by 
a factor of 11%. 



4. Results and discussion 

The significance map of the Cygnus region as observed by ARGO-YBJ using events 
with Npad > 20 is shown in Figure 1. For comparison, the 4 known TeV sources and 24 GeV 
sources in the second Fermi LAT catalog are marked in the figure. An excess is observed 
over a large part of the Cygnus region, which indicates a possible diffuse 7-ray e mission. A n 



analysis of the diffuse 7-ray emission using ARGO-YBJ data can be found in (IMa II2011I ). 
The highest significance value is 6.4 s.d. at (307.85°, 41.75°), consistent with the position of 
VHE sources MGRO J2031+41 and TeV J2032+4130. No evidence of an emission above 3 
s.d. is found at the location of MGRO J2019+37. 



-7- 



4.1. MGRO J2031+41 

The intrinsic extension of MGRO J2031+41 is determined by fitting the distribution 
of 9 2 for the events exceeding the background as shown in Figure 2, where 9 is the angular 
distance of each event to the position of TeV J2032+4130. Only events with N pa( i > 60 
are used in this fit, where N pat i is the number of fired pads. To fit the data, a set of 7- 
rays is generated taking into account the spectral energy distribution (SED), the intrinsic 
source extension, and the detector PSF. The extension is estimated, by minimizing the \ 2 
between the data and the generated events, between 0° and 1° with steps of 0.1°. Assuming 
the background spectral index —2.8, the intrinsic extension is determined to be a ex t = 
(O.2+Q2) . It is found that the dependence on the SED is negligible within the uncertainties. 
This result is consistent with the estimation by the MAGIC and HEGRA experiments, i.e., 
(0.083±0.030)° and (0.103±0.025)°, respectively. 

Assuming an intrinsic extension a ext = 0.1°, we estimate the spectrum o f MGRO 



J2031 +41 using the ARGO-YBJ data by a conventional fitting method described in ( IBartoli et al 



2011af ). In this procedure, the expectation function is generated by sampling events in the 
energy range from 10 GeV to 100 TeV and taking into account the detailed ARGO-YBJ 
detector response to the events assuming a power law with its spectral index as a parameter. 
We define four intervals with N pad of 60-99, 100-199, 200-499, and >500. The best fit to 
the SED is shown in Figure 3. The differential flux (TeV -1 cm -2 s" 1 ) in the energy range 
from 0.6 TeV to 7 TeV is 

" ?A 1.40 ± 0.34) x io-H(_|L)-2.83±o.37_ (!) 



dEdAdt v ' y lTeV 

The integral flux is 31% that of the Crab at energies above 1 TeV, which is higher than 
the flux of TeV J2032+4130 as determined by the HEGRA and MAGIC experiments, i.e. 
5% and 3% Crab u nit, respectively. Howe ver, this measurement is in agreement with the 



Milagro new result ( IBonamente et al.ll201ll ). also shown in Figure 3. 



The reason for the discrepancy between the fluxes measured by Cherenkov telescopes 
and extensive air shower arrays is still unclear. A contribution is expected from the diffuse 
7-ray flux produced by cosmic rays interacting with matter in the Galaxy plane. According 
to the mea surement of diffuse 7-ray flux from the Cygnus region using ARGO-YBJ data 



( )Ma Il201ll ). this contribution to the measured flux from MGRO J2031+41 at ene rgy above 1 



TeV is about 10%. Moreover, an estimate of the systematic error is described in ( lAielli et al. 



20101 ). With an incomplete list of the sources of systematics, such as the time resolution 
variation, event rate variation with environment parameters, and pointing error, the error is 
found to be less than 30%. Due to the limited angular resolution of the ARGO-YBJ detector, 
nearby sources could contribute to the flux. For TeV J2032+4130, this contribution must 
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be very small, because there is no source in the 3° x 3° field of v i ew as seen by HEG RA and 
MAGIC except TeV J2032+4130 itself fkharonian et all l2005[ [Albert et al] l2Q08h . Thus, 
the contribution from diffuse 7-ray emission, nearby sources, and systematic uncertainty are 
not enough to explain the discrepancy. 



4.2. MGRO J2019+37 



No excess above 3 s.d. is detected inside the Cyg OBI region, even considering extended 
sources. Taking into account the position uncertainty reported by Milagro, the bin with the 
maximum significance within 0.3° from MGRO J2019+37 is used to estimate the upper 
limits. The flux upp er limits at 9 0% confidence level (c.l.) are shown i n Figure 4 assumin g 
the SED reported in ( iSmith 1120101 ) and the extension a = 0.32° given in (lAbdo et al.ll2007bl ). 
respectively. At energies above 5 TeV, the ARGO-YBJ exposure is still insufficient to reach 
a firm conclusion. As regards the emission at lower en ergies, the A RGO-YBJ observation 
does not confirm the spectrum determined by Milagro ( iSmith II2010I ). Taking into account 
that the two observations differ in time by several years, this discrepancy might indicate a 
variation in the 7-ray flux of the source. 

The VERITAS experiment carried out a fine scanning of the Cygnus region. With a 
sensitivity of about 10% Crab unit, there was no significant signal found in the direction of 
MGRO J2019+37 (jWeinstein II2009I ). With a deeper survey correspondi ng to a sen sitivity 
of about 1% Crab unit, some faint sources were found in this region (lAliu 1 1201 ll ). The 
estimated flux is much weaker than that determined by the Milagro experiment. 

Considering the source extension a = 0.32° ±0.12° and the distance of the Cygnus region 
1—2 kpc, the source radius is estimated to be 4—15 pc, implying that the variation timescale 
should be longer than 13—49 years. The observation by the ARGO-YBJ experiment is about 
five years later than that by Milagro. A flux variation over the whole extended region cannot 
be completely excluded. If the flux variation were dominated by a smaller region in the source 
area, the picture could be more reasonable. In such a scenario, however, identifying MGRO 
J2019+37 as a PWN could be a dilemma because it otherwise should have a steady flux. 



5. Conclusions 

Since 2007 November the ARGO-YBJ experiment is monitoring with high duty cycle 
the northern sky at TeV photon energies. Using data up to 2011 August, we have observed 
the Cygnus region, inside which two bright VHE extended 7-ray sources have been detected 
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by the Milagro experiment. An excess with statistical significance of 6.4 s.d. is detected 
from the direction of MGRO J2031+41, consistent with the Milagro observation, but with a 
flux higher than that measured by HEGRA and MAGIC. The source location and extension 
are, however, consistent with those of TeV J2032+4130. It is not easy to assess the origin of 
this discrepancy. No signal from MGRO J2019+37 is detected, and the derived upper limits 
at 90% c.l. are lower than the Milagro flux at energies below 5 TeV. This result could be 
explained invoking a source variability, making difficult in that case to identify the source 
as a pulsar wind nebula. In conclusion, further observations and attention to the Cygnus 
region are needed since it is found to be complex in the VHE domain. 
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Fig. 1. — Significance map of the Cygnus Region as observed by the ARGO-YBJ experiment. 
The four known VHE 7-ray source are reported. The errors on the M GRO source positi ons 
are marked with crosses, while the circles indicate their intrinsic sizes fjAbdo et al.ll2007al lbl). 



The cross for VER J2019+407 indicates its extension (IWeinstein 1 12009). The source VER 



J2016+372 is marked with small circles without position errors. The small circle within the 
errors of MGROJ2031+41 indicates positio n and extension of the source TeV J2032+4130 
as estimated by the MAGIC collaboration ([Albert et al.ll2008l ). The open stars mark the 
location of the 24 GeV sources in the second Fermi LAT catalog. 
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Fig. 2. — Distribution of 8 2 for the number of excess events around TeV J2032+4130. The 
filled region is the best fit to simulated data. 
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Fig. 3.— Energy density spectrum from TeV J2032+4130/MGRO J2031+41 as measured 
by the ARGO-YBJ exp eriment (mage nta solid line). The spectral measurements of HEGRA 
(lAharonian et al.ll2005l ) and MAGIC ( lAlbert et al.ll2008l ) are also reported for comparison. 
The solid line and shaded area indicate the differential ene rgy spectrum and the 1 s.d. error 



region as recently determined by the Milagro experiment (jBonamente et al 



201 lh. The two 



triangles give the previous flux measurements by Milagro at 20 TeV (lAbdo et al.ll2007al ) and 
35 TeV JAbdo et al-lbopgh . 
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Fig. 4.— Upper limits to the flux from M GRO J2019+ 37 derived by the ARGO-YBJ ex- 
periment adopting th e spectrum given in (ISmith Il2010l ). The extension is assumed to be 
a = 0.32° as given in (lAbdo et al.ll2007bl ). The solid line and shaded area indicate the differ- 
ential energy spectrum and the 1 s.d. error region as determined by the Milagro experiment 



(lAbdo et al. 



(Smith 2010 ). The two triangles give the previous flux measurements by Milagro at 20 TeV 



2007al ) and 35 TeV flAbdo et all 120091 ). 



